Enhancement of variation of fundamental constants in ultracold atom and molecule 

systems near Feshbach resonances 
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Scattering length, which can be measured in Bose-Einstein condensate and Feshbach molecule 
experiments, is extremely sensitive to the variation of fundamental constants, in particular, the 
electron-to-proton mass ratio (m^/rrip or m^/ A.qcd, where Aqcd is the QCD scale). Based on 
single- and two-channel scattering model, we show how the variation of the mass ratio propagates to 
the scattering length. Our results suggest that variation oinie/mp on the level of 10~^^ ~ 10~^^ can 
be detected near a narrow magnetic or an optical Feshbach resonance by monitoring the scattering 
length on the 1% level. Derived formulae may also be used to estimate the isotopic shift of the 
scattering length. 
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Theories unifying gravity with other interactions sug- 
gest a possibility of temporal and spatial variation of 
the fundamental constants of nature (see e.g. review 
where the theoretical models and results of measure- 
ments are presented). There are hints for the variation 
of the fine structure constant, a — /Tic, strength con- 
stant of the strong interaction and fundamental masses 
in Big Bang nucleosynthesis, quasar absorption spectra 
and Oklo natural nuclear reactor data (see e.g. 12,13)0)13) 
. However, a majority of publications report only limits 
on possible variations (see e.g. recent review ^). 

The hypothetical unification of all interactions implies 
that variation of the electromagnetic interaction constant 
a should be accompanied by the variation of masses and 
the strong interaction constant. Specific predictions need 
a model. For example, the grand unification model dis- 
cussed in predicts that the quantum chromodynamic 
(QCD) scale Aqcd (defined as the position of the Landau 
pole in the logarithm for the running strong coupling con- 
stant) is modified as follows: 5Aqcd/^qcd ~ 34:Sa/a. 
The variation of quark and electron masses in this model 
is given by Sm/m ^ 706a/a. This gives an estimate for 
the variation of the dimensionlcss ratio 



(m/AQCD) a 



(1) 



The proton mass is proportional to Aqcd {mp ~ 
SAqcti)- Therefore, we have 

^(rrWrn^^g^fo (2) 
[me/mpj a 

This result is strongly model-dependent. However, the 
large coefficients in these expressions are generic for 
grand unification models, in which modifications come 
from high energy scales: they appear because the running 
strong coupling constant and Higgs constants (related to 



mass) run faster than a. This means that if these models 
are correct the variation of masses and strong interaction 
may be easier to detect than the variation of a. 

One can measure only variation of the dimensionlcss 
quantities, therefore we want to extract from the mea- 
surements variation of the dimensionlcss ratio rrie/mp. A 
number of limits on variation of rrie/mp and mq/ Aqcd 
( where is the light current quark mass with the de- 
pendence on the normalization point removed) have been 
obtained recently from consideration of Big Bang Nucle- 
osynthesis, quasar absorption spectra and Oklo natural 
nuclear reactor which was active about 1.8 billion years 
ago nil mm [Hill The present day limits have 
been obtained from the atomic clocks experiments. The 
most recent results have been presented in see also 
original measurements in [T^ . 

In this letter, we propose a new approach to measure 
the variation of electron-to- proton mass ratio /3 = rrie/mp 
with extremely high sensitivity. We show that observ- 
ables in atomic Bose-Einstein condensates (BECs) and 
Feshbach molecules experiments can be extremely sensi- 
tive to 13. This dependence can be understood as the 
molecular potential is essentially electronic, while the 
atomic mass is essentially baryonic. The high sensitivity 
results from the fact that the scattering length depends 
critically on the energy of the last bound state, while 
the variation of mass-ratio scales the overall vibrational 
energy, which lead to an stronger accumulative effect on 
the higher-lying bound states. 

Recent experiments have well characterized the atomic 
scattering length based on the mean-field shifts in 
BECs, coherent atom- molecule coupling and photo- 
association [T^, Feshbach fl^, and radio- frequency 
spectroscopy. In these works, scattering lengths are typ- 
ically determined to ~ 1% in the BEC-type experiments 
and to 10""* in the molecular spectroscopy experiments. 
In this paper, we theoretically evaluate the enhance- 
ment factor K, defined as the ratio of the fractional 
change of the scattering length to that of (3, namely. 
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5a/a = K{SP/f3). Our results show that huge enhance- 
ment factors oi K = 10^ ~ 10^^ can be obtained near 
narrow Feshbach resonances. This results can provide 
alternative stringent and independent tests on the varia- 
tion of fundamental constants. 

To quantify the dependence of scattering length on 
fundamental constants, we first consider a single-channel 
scattering model described by molecular potential V{r) 
which at large distances approaches the van der Waals 
interaction V{r) = —C^jr^. Here r is the atomic sep- 
aration and Cg is the van der Waals coefficient. This 
potential very well characterizes the interaction of ultra- 
cold neutral atoms. Atomic s-wave scattering length a in 
this potential is given by 



= a[l -tan(0-7r/8)], 



(3) 



where a = ci^ixC^jh})^!'^ is the mean scattering length, 
c = 2-3/2r(3/4)/r(5/4) w 0.47799 is a constant, is the 
reduced mass of the two interacting atoms and is the 
scattering phase shift. 

For most atoms, the associated molecular potential 
can support a large number of vibrational bound states 
and one can calculate the scattering phase shift semi- 
classically as |19i| 



h-^y/-2^V{r)dr, 



(4) 



where Vi is the classical inner turning point at zero energy. 
The number of vibrational bound state is = [(j)/7r+3/8] 
[Tgj i. where [x] is the integer part of x. 

Equation (O suggests that fundamental constants can 
influence the scattering length in two ways, either by 
varying the mean scattering length or by varying the 
scattering phase shift. Varying Eq. ||2J), we derive 



6a 6a {a — a)^ + 
a ~ 



a 



aa 



(5) 



The scattering length has dimension of length, there- 
fore, we have to specify units we will use. We are only 
interested in enhanced effects of the variation. In this 
case choice of units is not important since their variation 
can be neglected anyway. Note that the mean scattering 
length in Eq. Q can be presented as a = const. f3^^/'^aB 
where as is the Bohr radius. The phase (f) may be pre- 
sented as (j) ^ [An/PY''^ j yJ-V{r)/EHd[r/aB) where 
Eh — /meO^ is the Hartree energy. An is the number 
of nucleons in the nucleus. Therefore, it is convenient to 
use atomic units. We normalize all lengths to and 
energies to Eh- Performing variation of a and (j> we ob- 
tain that in atomic units they depend on variation of the 
electron-to-proton mass ratio P only: 



6a 



6(t> 



1 5/3 



(6) 



The variation of disa ppears since the classical mo- 
mentum p = y^—2iiV{r) is zero at the turning point. 
In the derivation of the above equations, we have used 
6{V/Eh) — 0. Indeed, the molecular potential is de- 
termined by the Coulomb interaction between electrons 
and nuclei. Therefore, in the non-relativistic limit all pa- 
rameters of this potential are proportional to the atomic 
units, there are no other natural units in this problem. 
The relativistic corrections bring some additional weak 
dependence on a (see e.g. which may be neglected 
here. Indeed, the relativistic corrections are of the order 
of Z^a^ and relatively small. For Cs Z^a'^ = 0.16, for Li 
Z^a^ = 0.0005 (here Z is the nuclear charge). 

Combining Eqs. lO and © and assuming the potential 
supports many bound states {4> « Nn 3> tt), we obtain 



6a _ Ntt [a- ay + a^ 6/3 
a 2 ad 



K 



6J3_ 
J' 



(7) 



where K is the desired enhancement factor. 

Note that this equation may also be used to estimate 
the isotopic shift of the scattering length if it is small 
{6P/P={An-A'j^)/An). 

Two enhancement factors are identified here: one from 
the number of bound states N and one from the resonant 
behavior of the scattering length ((a — d)^ + d?)/{ad). 
Note that we have large enhancement here in two limiting 
cases: very large a and very small a. 

The first enhancement factor is larger for heavier atoms 
with deeper molecular potentials as well as higher density 
of states. For example, singlet potential of Cs2 supports 
N « 150 vibrational states, while that of Li2 only has 
/V = 38 states. 

More remarkably, when the scattering length is large 
compared the mean scattering length, the second en- 
hancement factor approaches a/a, which can be very 
large near a potential resonance. In this limit, the bind- 
ing energy of the last molecular state is simply given by 
Eb — I2iic? (for a > 0). We can then write the varia- 
tion of the binding energy in the threshold regime as 



-Nn 



(a ~ a)2 6P 
ad (3 



(8) 



Atoms at low temperatures constitute an ideal sys- 
tem to measure the scattering length since the scattering 
length is the only parameter which characterizes the in- 
teraction properties. In the case of cesium, the scattering 
length in the singlet potential has been determined to be 
a — 280.37 a^ with a relative accuracy of ^ 10-4 ^ 
Based on the number of Cs2 vibrational states A^ « 150 
and the mean scattering length of a = 95 as, Eq. Q 
suggests an overall enhancement factor of AT ~ 400. 

Similar sensitivity oi 6a/ a — 2 x 10^'' on the singlet 
molecular potential of lithium-6 atoms is also obtained 
recently |l8j . Due to the relatively fewer molecular num- 
ber of A^ = 38 and the small singlet scattering length of 
a =45.17 as, the enhancement factor is only AT = 62. 
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For real atoms, the molecular potential is multi- 
channel in nature. This opens up new possibilities to 
modify the scattering lengths and to obtain higher en- 
hancement factors. In particular, when a bound state 
in a closed channel is tuned close to the scattering state 
in the open channel, Fcshbach coupling between the two 
channels can lead to a resonant enhancement of the scat- 
tering length according to 



a(l 



Eo — Er 



(9) 



where A is the scattering length near the resonance, a 
is the off-resonant scattering length, AE characterizes 
the Feshbach coupling strength, and Em is the energy of 
the bound state, which is supported by the closed chan- 
nel. We assume that the long-range asymptotic of the 
the closed channel potential is Vc — Ec ~ Cg /r^ and the 
asymptotic of open channel potential is Vo = Eo — Cq/t^. 
Note that the Feshbach resonance may have a finite decay 
width r„i, therefore, the resonance energy denominator 
is better described as Eo — Em + iVm/'^- However, the 
elastic width vanishes at zero kinetic energy of colliding 
particles and inelastic width is small. To simplify our 
consideration we will neglect Tm and other possible non- 
linear effects in the channel coupling. 

To determine the dependence of A on fundamental 
constants, we first assume that the molecular potentials 
of both channels and the Feshbach coupling are propor- 
tional to the atomic unit of energy, namely. 



,£'c Eo ,AE 

J^H J^H i^H 



(10) 



This assumption needs corrections if AE contains the hy- 
perfine interaction component ( see the variation of the 
hyperfine interaction in [l^'l. However, near the reso- 
nance the weak dependence of AE on fundamental con- 
stants is relatively insignificant and may be neglected. 
Variation of Eq. lO then gives 



5A 
~A 



Sa 
a 



A 



6E„ 



A Eo — E„ 



(11) 



We note that the energy of the bound state Em de- 
pends on the vibrational energy and therefore on the 
mass ratio. To estimate this dependence, we assume the 
molecular potential Vc{r) supports many bound states, 
and we may use a WKB approximation to estimate the 
molecular energy Em using equation 



Vc{r))dr ^ Mt: (12) 



where the integer M is the vibrational quantum number 
of the bound state and ri(ro) is the classical inner(outer) 
turning point for a total energy of Em- 



Equation (|12|l relates Em to the atomic mass m. Vary- 
ing Eq. H12I) , we can relate the level shift to the variation 
of the mass ratio f3. The result is 



Sf3 



^Em = -7, — 77r^-5- (13) 
2Tr p[Em) P 

where p{Em) = 1/-D is the density of states at the energy 
Em and can be estimated from the energy splitting D 
between adjacent vibrational levels. 

Combining Eqs. (|11I13() and using 0^/ = Mtt 3> tt, 
we show that near a Feshbach resonance, the change of 
scattering length is dominated by the resonance term and 
can be written as 



5A 
~A 



M {A 



1 



5(3 



Aa p{Em)AE P • ^^^^ 

There are three enhancement factors in Eq. 114|l : the 
first one is the vibrational quantum number of the closed 
channel bound state M; the second one depends on the 
scattering length and approaches A/ a near the resonance. 
These two factors resemble those from the single channel 
calculation, see Eq. j?)). With the Feshbach resonance, 
however, one is given the flexibility to adjust the second 
enhancement factor. The third factor of (/9(£'m)Ai<^)~'^ = 
D I AE in Eq. (14) is a new effect. It shows that higher 
energy resolution can be obtained by using a narrower 
Feshbach resonance and larger level shifts occur when 
the closed channel has a larger vibrational state energy 
splitting D. 

In the following, we take ultracold Li and Cs atoms as 
two examples to estimate the dependence of the molecu- 
lar binding energy or the scattering length on the mass 
ratio. For both atomic species, multiple Feshbach reso- 
nances exist and the vibrational energy splittings in the 
highest singlet channel near the lowest triplet scattering 
threshold aie D = h 10~20 GHz We wiU adopt 

p = D^^ = {h 10 GHz)^^ in the following estimation. 

For ^Li, a broad resonance at 834 G with a width of 
300 G (A£;=h 840 MHz) and a background scattering 
length of a =-1405 as 18] can be a candidate to reach 
very large scattering lengths. This resonance is induced 
by a closed channel bound state with M = 38 and molec- 
ular binding energy has been measured to ~ 1% near this 
resonance 18]. This measurement is performed deep in 
the quantum threshold regime with scattering lengths up 
to A ^ 1300 ttB- Here, we can rewrite Eq. H14I) in terms 
of the threshold binding energy Ei, = h^/2pA^ as 



5E, ^ _^^^iA_a 
Eb Aa 



1 



6(3 



(15) 



p{Em)AE (3 ■ 

Based on the Li parameters, we calculate the total en- 
hancement iohe K = —38 x (-4) x 12 = 1824, which is 
indeed much larger than the single channel result we ob- 
tained. The 1% binding energy measurement thus trans- 
lates into a mass ratio measurement with an accuracy of 
5 X 10-6. 
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Feshbach molecules of ^Li are also observed near a 
narrow resonance at 546 G with a width of 0.1 G 
(A£'=h 0.28 MHz) This resonance is induced by 

the same vibrational manifold with M — 38 and the 
background scattering length is a — 62 as- A measure- 
ment similar to |0| should be possible at ^ = 620 as, 
where we expect the enhancement factor to be as large 
asK = -38 X 8.1 x 36, 000 « -10^. 

An extreme case would be the narrow g-wsve Fesh- 
bach resonances observed in ^'^^Cs, which are induced 
by molecular states with four units of angular momen- 
tum. These resonances have a typ ical resonance width 
of ~ 5 mG {AE=h 3.5 kHz) ^ and can be induced 
by states in the closed singlet channel with M ^ 140. 
Since the threshold regimes of these resonances are ex- 
pected to be too narrow to access experimentally, we re- 
turn to Eq. (|14|l and calculate the variation of the scat- 
tering length. Assuming a scattering length measure- 
ment is performed at A — 1000 as near a narrow reso- 
nance and the background scattering length is 160 a^, 
we find the overall enhancement factor is as large as 
K = 70x (4.8) X (3 X 10^) « 10^. 

High resolution measurement near these narrow res- 
onances is, however, non-trivial since it requires excel- 
lent magnetic field stability. In the above example, a 
magnetic field precision of 1 mG is required to reach 
A = 1000 qb and a field instability and inhomogeneity 



of 0.01 mG are necessary to determine A to below 1%. 
Since an instability of 1 mG has been achieved near a nar- 
row resonance |2J| without magnetic field shielding, we 
are optimistic that 0.01 mG instability can be achieved. 
If so, a 1% measurement on the scattering length can be 
translated into a sensitive test of the mass ratio variation 
on the level of \Sf3/l3\ - 10"". 

Another possibility to obtain large enhancement fac- 
tors is by optical Feshbach resonances 25, 26]. In this 
case, the coupling strength can be externally controlled, 
and the bound state can be selected for higher sensitiv- 
ity on the mass ratio and very weak dependence on the 
magnetic field. Using the two-photon version of the op- 
tical Feshbach resonance for higher frequency resolution 
(100 Hz frequency resolution was reported in Ref. |l6j). 
we expect a coupling strength of AE =h 5 Hz can be 
achieved, which, in the case of cesium, can provide a large 
enhancement factor oi K ^ 10^^ and a test of /3 on the 
10~^^ level. Furthermore, the optical method allows one 
to probe molecular states with different dependence on 
(3. This can provide significant reduction of systematic 
effects. 

We thank Paul Julienne, David deMille and Jun Ye 
for stimulating discussions. V.F. is supported by the 
Australian Research Council. C.C. acknowledges sup- 
port from the Chicago Material Research Center. 



[1] J-P. Uzan, Rev. Mod. Phys. 75, 403 (2003); 

[2] M.T.Murphy, J.K. Webb, V.V. Flambaum. Mon. Not. R. 

Astron. See. 345, 609-638 (2003). 
[3] V.F. Dmitriev, V.V. Flambaum, J.K. Webb, Phys. Rev. 

D69, 063506 (2004). 
[4] S K Lamoreaux and J R Torgerson Phys. Rev. D 69 

121701(R) (2004). 
[5] A. Ivanchik, P. Petitjean, D. Aracil, R. Strianand, H. 

Chand, C. Ledoux, P. Boisse Astron. Astrophys. 440 , 

45 (2005). 

[6] S.G. Karshenboim, V.V. Flambaum, E. Peik. Atomic 
clocks and constraints on variations of fundamental con- 
stants. In "Handbook of atomic, molecular and optical 
physics" (Ed. G.W.F. D rake, Springer 2005), Chapter 
30, pp.455-463;"arXiv:physics/0410074 

[7] P. Langacker, G. Segre and M.J. Strassler, Phys. Lett. 
B 528, 121 (2002); See also X. Calmet and H. Fritzsch, 
Eur. Phys. J, 024,639 (2002), W.J. Marciano, Phys. Rev. 
Lett. 52, 489 (1984). 

[8] V.V. Flambaum, E.V. Shuryak, Phys. Rev. D65, 103503 

(2002) . 

[9] K.A. Olive et al, Phys. Rev. D66, 045022 (2002). 
[10] V.F. Dmitriev, V.V. Flambaum, Phys. Rev. D67,063513 

(2003) . 

[11] V.V. Flambaum, E.V. Shuryak, Phys. Rev. D67,083507 
(2003). 

[12] P. Tzanavaris, J.K. Webb, M.T. Murphy, V.V. Flam- 
baum, and S.J. Curran. Phys. Rev. Lett. 95 , 041301 , 
2005. 

[13] V.V. Flambaum, A.F. Tedesco, |nucl-th/060105'ol 



[14] H. Marion et al, Phys. Rev. Lett. 90, 150801 (2003). S. 
Bize et al., arXiv:physics/0502117 /I. J.D. Prestage, R.L. 
Tjoelker and L. Maleki, Phys. Rev. Lett. 74, 3511 (1995). 
S. Bize et al, Phys. Rev. Lett. 90, 150802 (2003). M. Fis- 
cher et al arXiv:physics/0311128 E. Peik, B. Lipphardt, 
H. Schnatz, T. Schneider, Chr. Tamm, S.G. Karshen- 
boim, arXiv:physics/0504101 

[15] N. R. Olaussen, S. J. J. M. F. Kokkelmans, S. T. Thomp- 
son, E. A. Donley, E. Hodby, and 0. E. Wieman Phys. 
Rev. A 67, 060701(R) (2003). 

[16] R. Wynar, R.S. Freeland, D.J. Han, 0. Ryu, D.J. Heinzen 
Science 287, 1016 (2000). 

[17] 0. Ohin, V. Vuletic, A.J. Kerman, S. Chu, E. Tiesinga, 
P.J. Leo, and O.J. Williams Phys. Rev. A 70, 032701 
(2004) 

[18] M. Bartenstein et al., Phys. Rev. Lett. 94, 103201 (2005). 
[19] G.F. Gribakin and V.V. Flambaum, Phys. Rev. A 48, 
546 (1993). 

[20] V. A. Dzuba, V.V. Flambaum, J. K. Webb , Phys. Rev. 

Lett. 82, 888, 1999; Phys. Rev. A59, 230, 1999. 
[21] P.S. Julienne, private communication. 
[22] K.E. Strecker, G.B. Partridge, and R.G. Hulet, Phys. 

Rev. Lett. 91, 080406 (2003) 
[23] J. Herbig, T. Kraemer, M. Mark, T. Weber, 0. Chin, 

H.-O. Nagerl, and R. Grimm, Science 301, 1510 (2003). 
[24] M. Mark, T. Kraemer, J. Herbig, 0. Chin, H.-O. Nagerl, 

and R. Grimm, Europhys. Lett. 69, 706 (2005) 
[25] P.O. Fedichev, Yu. Kagan, G.V. Shlyapnikov, and J.T.M. 

Walraven, Phys. Rev. Lett. 77, 2913 (1996) 
[26] M. Theis, G. Thalhammer, K. Winkler, M. Hellwig, G. 



5 



Ruff, R. Grimm, and J. Hecker Denschlag Phys. Rev. Grimm, and J. Hecker Denschlag Phys. Rev. Lett. 95, 

Lett. 93, 123001 (2004) 063202 (2005) 

[27] K. Winkler, G. Thalhammer, M. Theis, H. Ritsch, R. 



